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Abstract 

The service provided by mobile networks operated today is not adapted to spatio-temporal fluctua¬ 
tions in traffic demand, although such fluctuations offer opportunities for energy savings. In particular, 
significant gains in energy efficiency are realizable by disengaging temporarily redundant hardware 
components of base stations. We therefore propose a novel optimization framework that considers both 
the load-dependent energy radiated by the antennas and the remaining forms of energy needed for 
operating the base stations. The objective is to reduce the energy consumption of mobile networks, 
while ensuring that the data rate requirements of the users are met throughout the coverage area. 
Building upon sparse optimization techniques, we develop a majorization-minimization algorithm with 
the ability to identify energy-efficient network configurations. The iterative algorithm is load-aware, has 
low computational complexity, and can be implemented in an online fashion to exploit load fluctuations 
on a short time scale. Simulations show that the algorithm can find network configurations with the 
energy consumption similar to that obtained with global optimization tools, which cannot be applied 
to real large networks. Although we consider only one currently deployed cellular technology, the 
optimization framework is general, potentially applicable to a large class of access technologies. 
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I. Introduction 

Mobile networks are eonstantly growing, with the number of mobile subseribers worldwide 
exeeeding 7.1 billion in 2014 and the eurrent foreeast of 9.5 billion subseribers in 2020 Q. 
During this time, an even more drastie inerease is predieted for mobile data traffie, whieh is 
expeeted to grow by the faetor of 8. To keep paee with this trend, mobile network operators 
attempt to boost the network eapaeity by deploying additional antennas and base stations. This 
will lead to a dramatie inerease in the energy eonsumption, thereby eausing additional operational 
expenses (OPEX) and higher earbon dioxide (CO 2 ) emissions. Therefore, enhaneing the energy 
effieieney of mobile eommunieations systems has a positive double effeet: lower greenhouse gas 
emissions and redueed OPEX. 

Current mobile networks have been designed to provide the best possible serviee to the users 
at all times. Global network parameters and the network topology are largely statie, although, 
as pointed out in many studies (see for instanee [|3||, flUl, |[5l, dH), the traffie load fluetuates sig- 
nifieantly over time and spaee. Sueh spatio-temporal fluetuations ereate large eapaeity surpluses 
at times of low traffie demand, whieh in turn offers opportunities for energy savings through 
adaptation of the serviee supply to the aetual demand. However, in order to utilize the eapaeity 
surpluses for signifieant energy savings, it is essential to reduee the energy eonsumed by hardware 
and auxiliary equipment (e.g. eoolers), whieh is a dominant form of energy eonsumption in 
eurrent mobile networks. In faet, for a typieal network with today’s teehnology, base stations 
eonsume over 50% of the total network energy budget |[71. Most of the energy eonsumed by base 
stations is spent on powering hardware and auxiliary equipment. Erom this, we eonelude that 
signifieant energy savings ean be aehieved only by temporarily disengaging redundant hardware 
eomponents of a base station. Indeed, as pointed out by [l5]|, redueing the number of aetive base 
stations in periods of low traffie load offers a huge potential for energy savings. Switehing off 
base stations and introdueing a eooperation between different operators in an urban seenario is 
expeeted to lead to a reduetion in energy eonsumption of up to 29% (H. This effeet will beeome 
even more pronouneed in the future sinee, as aforementioned, the number of base stations is 
expeeted to inerease over the next years to meet the growing demand for wireless aeeess [|8l. 

A. Our contribution 

This paper deals with the problem of minimizing the overall energy eonsumption in the 
downlink ehannel of mobile (eellular) networks. By taking into aeeount the energy eonsumed 
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by hardware and auxiliary equipment, we address key shorteomings of most existing approaehes 
to the challenge of boosting energy efficiency of cellular networks. The underlying problem is 
of combinatorial nature because it essentially amounts to selecting a subset of network elements 
corresponding to the most energy-efficient network configuration, while providing the desired 
network coverage. More precisely, motivated by O, we formulate a combinatorial optimization 
problem to find a network configuration that consumes the least amount of energy, while 
satisfying traffic demands expressed in terms of minimum data rate requirements. In doing 
so, we balance different forms of energy consumption in an optimal manner by taking into 
account both the load-dependent energy used for transmission and the static energy consumed 
by hardware regardless of the actual load. Similar to IIH, the technology specific constraints are 
defined to capture the QoS requirements of the users. Although our optimization framework is 
generic in the sense that it can be applied to multi radio access technology (multi-RAT) systems 
by incorporating different RAT specific constraints, owing to the lack of space, our focus is on 
a single RAT according to the long term evolution (LTE) standard. 

The underlying combinatorial problem is in general hard to solve so, even for networks of 
moderate size, we follow a widely-known relaxation approach (see for instance ifT^ l and relax the 
problem by another problem as an intermediate step that can be solved efficiently. More precisely, 
we apply convex relaxation techniques to the constraints and approximate the objective function 
by a concave function. To try to minimize the concave function over the convex set, we use 
of a majorization-minimization (MM) technique [fTTll . Our algorithm exhibits major advantages 
over other existing schemes. In particular, the proposed algorithm is able to find good solutions 
in a relatively short time so that it can even be used to utilize fluctuations in traffic demand 
on a relatively short time scale. It can further cope with a variety of network elements such as 
cells, sectors or even antennas. As a result, we can take into account different sources of energy 
consumption to arrive at an energy efficient network configuration. In addition, the algorithm is 
able to balance the minimization of the static energy consumption, which requires the deactivation 
of as many network elements as possible, against the minimization of the load-dependent energy 
consumption, which calls for an appropriate load balancing to avoid highly loaded base stations. 
In order to ensure a broad range of applicability, the load-dependent energy consumption can 
be modeled as any concave or convex function of the load. 

We also discuss how our algorithm can be extended along two directions. First, we show 
how to incorporate techniques motivated by coordinated multi-point (CoMP) transmission llT^ . 
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which actually makes some relaxations superfluous and leads to further energy savings. Seeond, 
we use the load-based model of [[T3l for interferenee-eoupled eellular networks together with 
the framework of standard interferenee funetions lfT4l to show that larger energy savings are 
possible at the eost of slightly inereased eomplexity and overhead. 

B. Related Work 

Over reeent years, some researeh effort has been devoted to exploiting temporal and spatial 
redundaneies in wireless systems for energy savings. For instanee, Referenees ifTSl . llTbll . lUTl 
address the problem of finding an optimal number of base stations and eell site plaeements so as to 
minimize the overall energy eonsumption subjeet to some quality of serviee (QoS) requirements 
of the users. Assuming a wireless network based on time division multiple aeeess (TDMA), 
the objeetive of the study in lUSll is to minimize the overall expeeted energy eonsumption by 
optimizing the number of base stations and their loeations. The authors formulate the problem 
as a mixed integer programming problem and suggest using a simplex method together with 
the braneh and bound algorithm. The drawbaek of this approaeh is that, due to the TDMA 
assumption, the analysis does not earry over to systems with inter-eell interferenee, whieh is one 
of the major ehallenges faeed by designers of modem wireless eommunieation systems llT^ . 
m- Furthermore, braneh and bound methods may be slow EOl . whieh exeludes an applieation 
of these methods to real-time seenarios, even if the underlying problem is of moderate size. 

Referenees [flbl . If9l propose eentralized and deeentralized algorithms for wireless eommuni¬ 
eation networks to address the problem of base station seleetion in the presenee of traffle load 
fluetuations. Although the proposed approaeh seems to provide good solutions in reasonable 
time, it does not allow ineorporation of different sourees of energy eonsumption, whieh is of 
utmost importanee in modem networks eonsisting of hierarehieal struetures. In addition, the 
authors foeus on numerieal evaluations to justify the approaeh. No analytieal justifieation for the 
performanee of the proposed algorithms is given. 

The authors of ifTTl argue in favor of sleep mode teehniques eoupled with various network 
planning sehemes. A genetie algorithm is used to find energy-effieient network deployments, so 
the authors have developed a purely heuristie approaeh to put seleeted base stations into a sleep 
mode for energy-effieient network operation. In addition to the laek of any mathematieal justifi- 
eation, the main shorteoming of this work is that the proposed approaeh eannot ineorporate other 
radio teehnologies other than UMTS terrestrial radio aeeess network (UMTS: universal mobile 
teleeommunieations system). In eontrast, as mentioned before, our optimization framework is 
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general enough to be applied to multi-RAT seenarios, including the second, third and fourth 
generations of cellular networks [| 2 T]| . 

C. Notation and Paper Organization 

For a vector x G its ith component is Xi G M. Similarly, for a matrix X G its 

(«, j)-th component is Xij. Inequalities involving vectors, such as Xi > X 2 , are to be understood 
as component-wise inequalities. The set R+ denotes the set of non-negative real numbers, while 
R++ := R+ \{ 0 } is the set of positive real numbers. 

Given a matrix X G R^^^, we use x := vec(X) G R^^ to denote the vector obtained by 
stacking the columns of X. Note that the entries of x may be confined to take values on [0,1] 
or {0,1} depending on whether X G [0, 1 ]^^^ or X G {0, 

Definition 1 (fi-norm): For any vector x G R^ and matrix X G R^^^, their /o-norms |a;|o 
and |X|o are equal to the number of nonzero elements of x and X, respectively. For a scalar 
X G R, |x|o := 1 if X 7 ^ 0 and |x|o := 0 otherwiseQ 

The remainder of this paper is organized as follows. Section HI] introduces the underlying 
system model and in section |in| we outline the general problem to solve. In section |IV| the 
proposed algorithm to find solutions for our optimization problem is derived based on a worst- 
case inter-cell interference assumption. Section jV] presents how to explicitly take into account 
a more realistic inter-cell interference model. We present empirical evaluations of the proposed 
algorithm in section jVll 


II. System Model 

We consider the downlink channel of a multi-cell LTE network with an established network 
topology and a central network controller. The central network controller is responsible for 
collecting measurements, executing the proposed algorithm, and propagating updated network 
configuration parameters throughout the network. We assume that the network consists of L base 
stations. Each base station has multiple sectors (called cells in the following), and we denote 
the set of cells belonging to base station I by Si. The set of all base stations is denoted by C, 
and we use M. := Oi^c<Si to denote the set of all M cells in the network. The cell deployment 
is assumed to be dense enough so that coverage areas of different cells overlap. This implies 
that users can be served by different neighboring cells. 


'Although the Zo-norm is not a norm, we use the term “norm” as it is a common practice in literature. 
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A. Ensuring coverage via test points 

In order to ensure the desired eoverage anytime and everywhere in the eonsidered area, we 
impose eoverage eonstraints by adopting the eoneept of test points lf22l . 

Definition 2 (Test point): A test point (TP) is a eentroid of a pre-defined subarea that repre¬ 
sents an aggregated QoS requirement resulting from individual QoS demands of all potential 
users in this subareao Without loss of generality, we assume N TPs with the set of all TPs 
denoted by M := {1, 2, 

A eonsequenee of this definition is that small-seale fluetuations in QoS demand at the user 
level are averaged out at the TPs. These small-seale fluetuations must be eompensated by the 
lower layers of the protoeol staek (e.g. through adaptive modulation or eoding). 

Assumption 1: The QoS requirement for a TP eorresponds to the aggregated expeeted traffie 
over the respeetive area per unit time. This traffie requirement is expressed in terms of the 
minimum required data rate per TP. 

Assumption 2: If the minimum rate requirement of TP j is met, so are the requirements of 
the users in the assoeiated subareaE 

For serviees with no explieit data rate requirements (e.g. voiee ealls), we assume that they ean 
be supported if a minimum data rate per serviee request is ensured. By Assumption [U eaeh 
TV j E fif is assigned rate requirement Vj, and we eolleet the rate requirements of all TPs in 
the veetor r = [ri, r 2 ,..., r^v] G ^++- In general, a TP ean be assigned to any eell, and an 
assignment should be understood as follows. If TP j G M is assigned to eell i G Ai, then all 
users in the respeetive subarea assoeiated with TP j are served by eell j. The assignment of 
the TPs to the eells is subjeet to optimization in this paper. We use X = [xifi G {0,1}^^^ to 
denote the assignment matrix where Xij = 1 if TP j is assigned to eell i and Xij = 0 otherwise. 


Assumption 3: While eaeh TP is assigned to exaetly one eell, eaeh eell ean serve multiple 
TPs, and the set of TPs served by eell i under assignment X is denoted by Affix.) c Af. 

We point out that this assumption has been widely used in previous studies ^22} . [|^ , ifT^ . and it 
is valid throughout the paper exeept for Seetion llV-Ci where it is shown how to inelude seenarios 


test point becomes a user if it represents a QoS requirement of one particular user, in which case the subarea is a point 
corresponding to the position of this user. 

^The smaller the area represented by each TP, the better is this approximation. However, smaller areas imply an increased 
number of TPs, and the computational complexity of the proposed algorithm grows. 
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in which each TP can be served by multiple eells. Note that if A/i(X) = 0 for some i G A4, 
then eell i ean be deaetivated for energy savings beeause no TP is assigned to eell i. In eontrast, 
if Afi(X) 7^ 0, then eell i is aetive, and eaeh TP eonneeted to it induees some amount of cell 
load. 

Definition 3 (Cell Load): Given the assignment x := vec(X), the load of eell i, denoted by 
Pi{x) G [0,1] or simply pi for notational simplieity, is defined to be the ratio of the number of 
resouree bloeks requested by TPs served by eell i E M. io the total number of resouree bloeks 
Bi available at this eell0 

We use p -.= [pi,..., pmY ^ [0^ 1]^ 1° denote the veetor of all eell loads. From the definition 
of eell load, we have the following: 

Fact 1: The load at eell i satisfies pj > 0 if and only if (iff) eell i serves at least one TP 
B. Spectral efficiency and resource usage 

The optimal assignment of TPs to eells is strongly influeneed by the speetral effieieney of the 
eorresponding links. For the analysis in this paper, we adopt an OFDMA-based model for the 
speetral effieieney that is widely used in the literature iTTSlI . Il2^ . [1241 . The speetral effieieney also 
depends on radio propagation properties. Therefore, we assoeiate to eaeh TP a path-loss veetor 
and write the path-loss veetors of all TPs as eolumns of the path-loss matrix G = [pifi G 
where p* j eaptures the long-term path loss and shadowing effeets for a radio link from eell i to 
TPj. 

Assumption 4 (Reliable path-loss estimates): A reliable estimate of G is available at the een- 
tral network eontroller. 

Remark 1: The problem of reliable estimation and traeking of the path-loss matrix is out of 
the seope of the paper. However, the matrix eaptures only long-term fading effeets, so reliable 
estimates of G ean be obtained and traeked in praetiee. Promising algorithmie solutions to this 
estimation problem are for instanee presented in [l25l . Moreover, in network planning problems, 
knowledge of G is a very eommon assumption in the literature [fT3l . ifT^ . [l22l . 

Now we are in a position to define the signal-to-interferenee-noise-ratio (SINK) 

M+ between eell i E M. and TP j G AA by lfT3l . ifT^ : 

7ij(p) = (1) 

^keM\{i} 9k,j Pk + 

''Note that Bi can also be interpreted as the total bandwidth available at cell i, in which case pi is expressed in terms of the 
fraction of required and available bandwidth. 



where Pj > 0 is the transmit power per resource block of cell i and cr^ > 0 is the noise power per 
resource block. Accordingly, the link spectral efficiency ojij : M:||^ —>■ R+ (in bits per resource 


bloc la for the link from cell i to TP j is given by lf2^ 



( 2 ) 


where and G M++ are suitably chosen constants, referred to as bandwidth 

and SINK efficiency, respectively. These constants depend on the overall system design, which 
includes the choice of scheduling protocols and multi-antenna techniques. The choice of these 
constants has no impact on our results, so they are assumed to be arbitrary and fixed throughout 
the paper. For realistic values of these constants, we refer the interested reader to IfT^ . If2^ . 

From (l2l), we can easily see that the necessary number of resource blocks bij at cell i to serve 
TP j with data rate Vj is equal to bij = > 0. In addition, following Definition [3l the load 

at cells can be computed by the following system of non-linear equations 



(3) 


Remark 2: In practice, cells need to reserve some fraction of their resource blocks for signal¬ 
ing. If cell i has B* resource blocks in total, and it needs to reserve Oi > 0 of its resource blocks 
for signaling, then the resource blocks at cell i available for allocation to TPs are Pj = B* — ai. 

For a fixed assignment X cell load p in (jS]) can be efficiently computed by means of fixed- 
point algorithms (c.f. Sec. However, the assignment of TPs to cells is the main subject of our 
optimization problem and thus we cannot evaluate dH) easily. In order to keep the complexity of 
the optimization problem tractable, we lower bound the spectral efficiency. 

Assumption 5 (Worst-Case Interference): We have the worst-case interference scenario if all 
cells are fully loaded, i.e. p = 1. 

Unless otherwise stated, we assume the worst-case interference, which results in a lower bound 
on the link spectral efficiency u:i^j{p) > Uij := for every p e [0,1]^. In general, this 

bound diminishes gains in energy savings when taking into account the energy consumption 
of hardware, and we show in Section |V] how to incorporate the actual link spectral efficiency 
to improve the energy savings. Nevertheless, having fully loaded cells as in Assumption |5] is 

resource block is defined as a portion of the available time-frequency plane spanning a number of consecutive OFDM 
symbols in the time domain over a number of sub-carriers in the frequency domain. 





desirable beeause it has been proven in that full load (i.e. p = 1) is optimal with respeet 
to the transmit energy eonsumption (see also [l27ll l. 

Remark 3: The worst-ease interferenee assumption eannot exploit the full potential for energy 
savings, but the assumption is of high praetieal relevanee beeause it is an effective way to avoid 
coverage holes as a result of deactivating cells based, for instance, on imperfect information. 

C. Energy consumption model 

In contrast to most work in literature, we consider a model for the energy consumption of a 
base station and its cells that takes into account not only the cell load-dependent transmit energy 
radiated by antennas, but also the remaining sources of energy consumption that are independent 
of the cell load as long as the cell/base station is active. 

Definition 4 (Active base station/cell): Consider a particular base station I G C and its cells 
i E Si. Let Pi E [0,1] be the load of cell i. We say that a cell i is active iff pi > 0 and that base 
station I is active iff one of its cells is active, i.e. Pi > 0. If a cell or base station is not 

active it is said to be inactive. 

With Definition |4] we are in the position to define the energy consumption of a base station. 

Definition 5 (Energy consumption): Given a TP assignment X inducing a cell load p, the 
energy consumption Ei(p) > 0 of base station I is defined to be the power that the respective 
base station consumes per unit of time, where Ei(p) = 0 iff base station I is inactive. 

The function Efip) depends on the hardware setup of the base station, but it can be split into 
three parts in general: 

(i) The static energy consumption of the base station Cl >0 (due to shared hardware between 
sectors e.g. cooling, power supply, etc.), 

(ii) The static energy consumption Ci > 0 (i E Sfi of its active cells (e.g. due to power 
amplifiers, signal processing units, etc.), and 

(iii) the load-dependent dynamic energy consumption of its active cells ffipfi (i E Si), where 
fi : [0,1] —>■ M+ is a given continuous function relating the energy consumption to the 
corresponding cell load. 

By these definitions and Fact [H Efip) is a discontinuous function of the cell load, and we 


have 




cells i E Si serve no TP 
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where iS^^active C Si is the set of of active cells of base station 1. Therefore, the total energy 
consumption in a network, which is the accumulated energy consumption of all active base 
stations, yields 

= = Y1 + + )• 

lec i&c V i£Si / 

For concreteness, we make the following assumption throughout the paper (see also Remark Hj) 
Assumption 6 (Concave dynamic energy consumption): fi : [0,1] —>■ M+ {i G A4), is concave 
and continuously differentiable. 

In particular, this assumption is satisfied by a linear dependency of the base station energy 
consumption and the cell load reported in current studies such as [[28ll . 

Remark 4: In fact, the load-dependent dynamic energy consumption can also be assumed to 
be a convex function of the load. Moreover, we could even assume that it is a sum of convex and 
concave functions. The optimization framework presented in this paper can be straightforwardly 
extended to cover these cases. 


III. Problem Statement 

Spatio-temporal redundancies in coverage and capacity resulting from day-time fluctuations 
in traffic demand present great opportunities for energy savings by deactivating redundant cells 
at times of relatively low traffic demand. Indeed, if the traffic demand decreases, some or all 
entries of the rate requirement vector r G become relatively small, which can be utilized to 
reduce the total energy consumption by minimizing the cost function in dH) subject to different 
constraints that follow from the system model and ([3]). Formally, the problem under consideration 
can be stated as follows (note that the complete set of equations is referred to as ([5])): 


min. d 


+(ci Pi 0+ fiipi)) 1 


(5a) 

l£C \ 


0 / 



r,- 

S- t- R ~ = Pi 

i G M 

(5b) 





^^4 

= 1 


j eAf 

(5c) 

ieM 





pi £ [0, 1] 


i G M 

(5d) 

Xij G {0, 1} 


i e M,j e Af, 

(5e) 
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where the optimization variables are Xij and pi {i E Ai, j G M). In partieular, Assumption [3] is 
captured by (l5^ together with (l5el) . Constraints (l5bl) and (l5dl) . in contrast, ensure that the cell 
load is in accordance with Definition [3l 

To ensure feasibility of the above problem and to show the effectiveness of our approach, 
we consider scenarios where the rate requirements of TPs are sufficiently low for a reasonable 
amount of redundancies that allow for deactivation of cells. Moreover, if the traffic requirements 
in the system are sufficiently low or the number of cells is sufficiently large, p* is expected to 
be sparse with zero entries specifying cells that can be deactivated. 


IV. Energy-Efficiency Optimization 


The difficulty of problem Q lies in its combinatorial nature. In fact, it can be shown that the 
problem is related to the classical bin-packing problem, which is known to be NP-hard [|29l . 
Consequently, the complexity is expected to grow exponentially with the number of cells. On 
the positive side, problem ([5]) has a special structure that can be exploited by majorization- 
minimization techniques [fTTll . which have been widely used in recent years to tackle various 
problems in compressed sensing [[30ll and machine learning ifSTI . 

Instead of finding a global solution to dS]), we will pursue a less ambitious goal. We apply 
the majorization-minimization techniques mentioned above to develop a low-complexity anytime 
algorithm that has a strong analytical justification. This algorithm is expected to provide good 
results (in terms of low energy consumption) with low-complexity. To this end, we reformulate 
problem Q to pose it in a more tractable form. Eirst, we observe that each load pi is, in fact, 
a function of X (c.f. Definition [3] and (l5bl) l. We can therefore modify the problem to have 
only X as an optimization variable. Recall that, if at least one TP is served by cell i ( i.e., 
— 1)’ then it follows from Pact [H that the cell load at cell i is non-zero and |pj|o = 1. 
Hence, the objective function in (l5al) can be equivalently written as 


I 

i&c 

= E 

lec 


a Pi 

ieSi 


+ ^2 + Mpi)) 

Q ^^* 5 / 


Cl 


ieSi j&M 


+ ^ I ei 

Q 






Y, ('"'It 

lec \ 


iGSi 


( 6 ) 
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where s* := vec(Sj) with Sj G {0,1}^^^ being a matrix of zeros, exeept for its ith row, whieh 
is a row of ones. Similarly, ti := vec(Ti) with G {0,1}^^'^ is a matrix of zeros, except for 
its rows i & Si, which are rows of ones. 

Definition 6: Given the assignment x and the load dependent energy consumption ffipfix)) 
of cell i with Pi{x) = dSH)), we define the function fi : [0,1]^^ —)■ M_,_ ; 

ME 

Considering Definition!^ and using Pi < 1 (see Definition |3]) in (l5bl) . we arrive at an equivalent 
problem given by 


min. ( Cl \tjx\^ + y] (a sf* ^ + Mx)] 

) 

(7a) 

lec \ ieSi ' ' / 



r,- 

s. t.: > — XiA <1 

j&N 

i G M 

(7b) 

^ 

j efif 

(7c) 

i&M 



Xij G {0, 1} 

i e M,j e Af, 

(7d) 


where the assignment variables Xij (i e fiA,j e J\f) are the only optimization variables. 


A. Problem relaxation 

To obtain an optimization problem that is computationally tractable, we first relax the binary 
constraint (l7dl) tqj 


Xi^j G [0,1] , Vi G Ad, Vj G M. 


( 8 ) 


The above makes all constraints convex, so now the only problem is the objective function, 
which is not continuous due to the Zo-norm. We also note that by Assumption 0 and Definition 
the load-dependent term ffix) in the objective function dTai) is concave and continuously 
differentiable for a; G [0, 1]^^ since these properties are preserved under a composition with a 
linear function [[32ll . If^ . To address the non-continuity of the Zo-norm, we consider the following 
relation [|^ : 


l^lo 


log(l + e-i) 


(9) 


®This relaxation together with l7^ leads to a communication scenario where multiple cells serve one TP. A more detailed 
discussion on the implications is presented in Sec. lIV-Cl 
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By using ([91) and the non-negativity of Si, ti, x, the cost function in (TTal) can be equivalently 


written as 




lec 


ieSi 


lim > 

£->■0 ^ 


l€C 


log 1 + e tfx) ^ 
T^ + Z^ 

ieSi 


log(l -f e 




( 10 ) 


log(l -f e 


We can therefore obtain an approximation to problem ([5]) by replacing the objective function by 


the right-hand side of (fT^ for a sufficiently small but fixed e > 0. More precisely, for some e > 0, 
the objective is to find a matrix X G [0,1]^^^ or, equivalently, a vector x = vec(X) G [0,1]^^ 


that solves the following problem 

log(l + e-^ 


mm 


log(l -f e” 



E 

ieSi 


log(l -f e 


^ sTx) 


log(l -f e 


-n 


+ fi{x) 




5 ^ Xi,3 = 

i&M 


(11a) 
i G M (11b) 
j G AA (lie) 


Xij G [0,1] i e MJ e M. (lid) 

Solving problem (fTTI) is not straightforward because we need to minimize a non-convex function 
over a convex set. Fortunately, Reference [f30l presents an optimization framework based on the 
majorization-minimization (MM) algorithm ifTTll to handle problems of this type. The framework 
can be used to decrease the value of the objective function in a computationally efficient way. 
For completeness, we the reader can find some details of the MM algorithms in the appendix. 


B. Majorization-minimization (MM) algorithm 

For notational convenience, we define q := iog(i+e-i) ’ these 

definitions in (II lal) to simplify the objective function (ignoring unnecessary constants): 

h: X a; H- h{x) = ^ (q log(e -f tfx)) + ^ log(e + sfx) + fi{x)'^ , (12) 

l€C ieM 

where X C R^^ is the closed convex set of points satisfying the constraints (II lbl) - (ll Idl) and 
we have used the fact that A4 = U;g£iS;. Since /* is concave and continuously differentiable by 
Assumption!^ so is the function in ([T^ for any e > 0. Therefore, according to the explanations 
in the appendix, we can use the following function 

g : X X X ^R: {x,y) ^ h(y) + Vh{y)^{x - y) 
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as a majorizing function of (fT^ . where the gradient ean be easily ealeulated: 




e* 


e + sfx 


+ v/.(i; 


( 13 ) 


lec ~ ' ‘ ieM 

Thus, updates of the MM algorithm take the form (see the appendix) 

;("•+!) g argmin5f(a;, 




xex' 


tjx 


= arg min > Ci ^ , 


E 


^ i£M 

(0) 


T ~ 

Sf X 


e* 


e + 




X 


(14) 


for some feasible starting poinj^ G X. In words, the MM algorithm solves iteratively a 
sequenee of eonvex optimization problems. For the ehosen majorizing funetion, the problem to 
be solved in every iteration is a linear programming problem (LP), whieh ean be typieally solved 
effieiently with standard optimization tools. 

As diseussed in the appendix, the sequenee C X for some x^^'^ G X generated 

by (fT4l) produees a non-inereasing sequenee of objeetive values. Therefore, as 

n —)■ oo, we expeet the eorresponding sequenee of assignment matriees (note that 

=; vee(X("))) to evolve towards network eonfigurations with low energy eonsumption. 

We stop the algorithm if the improvements in the objeetive value are small enough in the 
sense that for some suffieiently small e* > 0, the following eondition is met 


h{x^^'>) - h{x^^+^^) < e\ (15) 

Upon termination, the resulting assignment matrix G [0,1]^^^ needs to be mapped to a 
matrix X* G {0, in order to obtain a feasible point to the problem in ([5]). For this purpose, 

we use the heuristie deseribed in Alg. [IJ The main idea is as follows. We start by rounding the 
entries x^J to the elosest integer, and then we eheek if the obtained assignment matrix is part 
of the set X. Otherwise, we aetivate additional eells and eonneet TPs to them. By using the 
standard LP solver of CPLEX, in our simulations most entries of the matrix G [0, 
are typieally either zero or one, so the rounding operation rarely results in a violation of a 
eonstraint (but we emphasize that this is not guaranteed to be true in general). 

For eonvenienee, we summarize the eomplete approaeh in Alg. [21 


^In our experience a good starting point is derived from a feasible assignment matrix obtained by connecting each TP to the 
cell providing the strongest received signal strength. 
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Algorithm 1 Heuristic to map [0, jo, 

Input: M, M, set of constraints X 2 representing (lllbl) and (l5el) 

Output: final assignment matrix X* 

1 : initialize: set of assigned TPs A = % and final assignment matrix X* = 0. 

2 : for all i G Ad, j G A/" do 

3 : if x-j G {1} then 

4: x*j = x-j and A = AVJ{j]. 

5: end if 

6 : end for 

7: Define set B = G (0,1) |V1 G M,Vj G A/'\Al|. 

8: while B ^ ^ do 
9: (l,j) = argmax^jjS} 

10 : if x*j := 1 -> X* G ^2 then 
11 : x*j = 1 and A = AU {j}. 

12 : B = B\{xi"^ \Wi e M} 

13: else 

I* i3 = B\{i.'”'} 

15: end if 

16 : end while 

17: for all j ^ Al do 

18 : activate closest non-active cell i which yields x*j := 1 —X* G ^2 and assign x*j = 1. 

19: A = AU{j}. 

20 : end for 


C. Serving a test point with multiple cells 

By Assumption [3l each TP is restricted to be served by exactly one cell. This strict limitation 
introduces the non-convex constraint (l5el) to the optimization problem in Q, which motivates the 
relaxation ([8]) and the heuristic mapping introduced in Alg. [U To avoid these heuristic approaches 
for which we are not guaranteed to find solutions, we assume in this section that each TP can 
be served by multiple cells. This assumption is implemented by using ([8]) directly instead of 
(l5el) . As a result, there is no need for any relaxations of the constraints or the use of heuristic 
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Algorithm 2 Network reconfiguration for improved energy efficient operation 

Input: set of TPs, set of cells, constraints 

Output: optimized network configuration according to X*. 

1 : initialize with a feasible point. 

2 : repeat 

3: compute by solving (fT4l) 

4: increment n 

5: until (fTSl) is valid 

6: use Alg. [Uto map to X* G {0, 

7: connect the TPs to cells according to X*. 

8: deactivate all cells no TP is connected to. 


mappings such as that in Alg. [TJ We only need to approximate the cost function as done in 
(II lal) and apply the MM algorithm to the resulting optimization problem, and we note that these 
operations have a strong analytical justification. 

The assumption of multiple cells serving one TP has a practical interpretation when considering 
Definition [21 It means that cells can serve only a fraction of the traffic generated in the area 
corresponding to some TP. In other words, we do not use a all-or-nothing approach, where cells 
should serve either all users or no users in the area corresponding to a TP. 

V. Load-Aware Energy-Efficiency Optimization 

The model presented in Sec. jn] assumes the worst-case interference in a fully loaded system, 
which leads to a lower bound on the link spectral efficiency (c.f. Assumption [5]). As pointed 
out in Remark [3l the main rationale behind this approach is the need for avoiding coverage 
holes when network elements are deactivated. The price is a sub-optimal performance in terms 
of energy efficiency because the interference is overestimated, and therefore users may use more 
resource blocks than required to keep their minimum data rate requirements. An immediate 
consequence of this is that more cells are activated than are necessary for meeting the minimum 
rate requirements at the TPs. In this section, we extend the optimization problem in (fTTI) to 
incorporate more precise estimates of the load induced by a given user-cell assignment, which 
is not a trivial task because it involves load computation (with fixed assignments) that requires 
the solution of a system of nonlinear equations ifT^ . |[T3ll . [l34l (note that we can easily estimate 
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the link spectral efficiency from the load by using ([21)). 

In what follows, we propose an approach that typically yields good approximations of the 
true link spectral efficiencies. The idea is to use a two-step alternating iterative scheme: 

Step 1 Compute the link spectral efficiency ^i,j{p) defined in ([21) for the load value 

obtained in the previous iteration of Step 2 of the algorithm (in the first iteration of the 
algorithm, we can use the worst-case spectral efficiency) and solve Problem (fTTI) with 
these (fixed) link spectral efficiencies to obtain an TP-cell assignment X. 

Step 2 For the TP-cell assignment obtained in Step 1, compute the load induced by this assign¬ 
ment. 

Regarding the load computation in Step 2, we use the fact that the load p induced by a given 
assignment X is a fixed point of the following standard interference mapping (see llT4l . [[^ and 
the references therein for further details): 

^ : p ^ [h{p) • • • Im{p)Y , 


where 


Ii{p) := min ^ 


^i,j ^i,j 


logs (l + 


r is a large constant and Ajj := • Since JT is a standard interference mapping and /j(p) is 

bounded above, we conclude that the fixed-point always exists and is unique ffT4ll . [l36l . Moreover, 
efficient iterative methods are known to approach the fixed point with an arbitrary precision [[T4l . 
061 . We summarize the heuristic proposed in this section in Alg. [3l 


VI. Numerical Evaluation 

In the following we present a numerical evaluation of the performance of the proposed 
algorithm in different networks. We start by outlining the basic simulation scenario followed by 
a comparison with two reference schemes with respect to the energy savings and computational 
time. Next, we present the ability of the proposed algorithm to incorporate a variety of different 
base station energy consumption models. Finally, we show the performance gains achieved by 
applying Alg. [3] from Sec. |Vl 
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Algorithm 3 Load-aware energy minimization 

Input: Worst-case spectral efficiency = a;(l). Maximum number of iterations Z. 
Output: Network configuration with low energy consumption. 

1 : for n = 0 : Z do 

2 : Use to construct Problem ([5]). 

3: Use Alg. [2] to obtain and remove deactivated cells from the set of cells to be 

considered in subsequent iterations. 

4: Compute the new link spectral efficiency for the assignment X^"^) by computing the 

fixed point of the standard interference mapping J. 

5: end for 

6: Return the network configuration resulting from X^^^. 


A. Basic Simulation Scenario 

The simulated network is located in a square-shaped area of size 2km x 2km, where L base 
stations are placed at locations chosen uniformly at random. Unless stated otherwise, each base 
station has three cells directed at 0°, 120° and 240° respectively. Traffic generated by users is 
represented by N TPs on an irregular grid. Hence, each TP represents the traffic requirements of 
an area of different size. To obtain spatially varying traffic requirements, we use the following 
traffic model in each run of the simulations. We define three circular hot-spot areas with centers 
chosen uniformly at random within the area. There are two types of TPs: “hot-spot TPs (HTP)” 
and “standard TPs (STP)”. Each TP in the simulation has probability 0.3 of being a HTP and 
probability 0.7 of being a STP. While the position of STP is chosen uniformly at random within 
the whole area, a HTP can be assigned uniformly at random to one of three hot-spot area. Its final 
position is determined in polar coordinates by sampling the distance from the hot-spot center 
from a normal distribution and the angle from a uniform distribution. We use a wrap around 
model to avoid boundary effects and determine the location of TPs to be placed outside the 
square-shaped area. The data rate requirements of TPs are derived from a normal distribution 
with /id = 128 kbps and variance aj = 32 kbps^ with a lower bound of 1 kbps. The signal 
attenuation for links between cells and TPs follows the ITU propagation model for urban macro 
cell environments with a horizontal antenna pattern for 3-sector cell sites with fixed antenna 
patterns lf37l . 
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Unless otherwise stated, we use the following simulation parameters: e* = 10“^, e = 10“^, 
Bi = 20MHz, Pi = 4QdB, t^sinr = 1, ^bw = 0.83, Cj = 500W and = 280W. The values of 
the last six parameters have been ehosen to mimie the behavior of eommereial LTE systems. 
Furthermore, we use fi{pi) = 564 pi to model the load-dependent energy eonsumption, whieh 
is a value similar to the dynamie energy eonsumption of current macro cells with 6 transmit 
antennas [[2^ . 

The proposed algorithms are compared with a solution of the original problem in ([5]) and, where 
possible, with the centralized cell zooming approach from ifT^ . The solution to the problem in 
([5]) is obtained by using Matlab 2013a in combination with IBM’s CPLEX on a Intel Core i7 
PC with four cores. As shown later in this section, the computational time to solve Q grows 
fast with the problem size. Therefore, to solve the problem in dS]) in a reasonable time for 
comparison purposes, we confine our attention to small networks with M = 102 cells (L = 34 
base stations) and N = 100 TPs, unless otherwise stated. We obtained the 95% confidence 
intervals depicted in the figures by applying the bias corrected and accelerated bootstrap method 
ll38l to the outcome of 100 independent runs of the simulations. Results related to the overall 
network energy consumption will be normalized to the energy consumption of the network when 
all cells are active and fully loaded. 

Definition 7 (Normalized network energy consumption): Given a TP assignment X inducing 
cell load p and given the resulting network energy consumption E(p), the normalized network 
energy consumption is defined to be 

E(p) E(p) 


Ap) ■■= 


-^(1) ZlzeC+/i(l)) 


where the term in the denominator is the energy consumption for a fully loaded system (p = 1). 

We refer to the sparsity supporting majorization-minimization algorithm as “sMM” and to 
any algorithm that solves Q directly as “MIP” algorithm (MIP: mixed-integer programming). 
We refer to solutions obtained by the centralized cell zooming algorithm in lIT^ as “cCZ”. The 
alternating approach proposed in Sec. |V]is referred to as “alternating sMM” algorithm. 


B. Computational performance comparison between sMM, cCZ and MIP 

The cCZ has limited capability to incorporate different energy consumption models and base 
stations with several sectors, so we confine ourselves to a simple base station model. We assume 
a homogeneous network model under which all base stations have only one omni-directional 




20 



0 ^-^-^^^^^^^^^ 

100 200 300 400 500 600 700 800 900 1000 


number of TPs 

Fig. 1: Comparison of normalized network energy consumption obtained with the sMM 
algorithm, the cCZ algorithm and the solution of the MIP problem for increasing number of 
TPs. Normalization with respect to the network energy consumption for a fully loaded system 
(p = 1) when all cells are active. 


cell, and all base stations have the same energy consumption model. More precise, we use 
\C\ = M = 100, |iS;| = 1 and dH) with q = 500, e, = 280, fi{pi) = 0 (/ G £, z G M). 

To show trends, we start with the standard setup described above, and we gradually increase 
the number of TPs in the system. Fig. [U shows the normalized network energy consumption. As 
expected, the normalized network energy consumption for all three algorithms increase as the 
number of TPs increases. This is intuitive because additional TPs add extra rate requirements 
that increase the total system load, which in turn reduces the redundancy in the network to be 
exploited for energy savings. The proposed sMM algorithm as well as the MIP algorithm provide 
network configurations that exhibit much smaller normalized network energy consumption when 
compared with the network configurations obtained with the cCZ algorithm. The smallest energy 
consumptions are achieved with the MIP algorithm, which outperforms the proposed sMM 
algorithm. For the scenario with 200 TPs, the sMM algorithm results in normalized network 
energy consumption of 12% on average. For the same number of TPs, the average normali z ed 
energy consumption under the cCZ and MIP algorithm are 49% and 7%, respectively. Similarly, 
for 1000 TPs, the resulting average normalized network energy consumption of 31% for the sMM 
algorithm is still larger than the 21% normalized energy consumption corresponding to the MIP 
solutions. However, it is still much smaller than cCZ with 88% normalized energy consumption. 
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These results emphasize that the sMM algorithm is a suboptimal heuristic, which is able to 
find network configurations consuming low energy. Even though the resulting network energy 
consumption is not globally optimal, it shows much larger energy savings than the comparison 
scheme cCZ. 

The main advantage of the proposed sMM algorithm is its fairly low computational complexity, 
which is directly affecting the time required to obtain an optimization result. Fig. [2] depicts the 
normalized time needed to obtain the results of Fig. [IJ This time is normalized with respect to the 
computation time of the MIP algorithm with 100 cells and 100 TPs. The sMM algorithm always 
provides results in a substantially shorter time than the MIP algorithm. Even for a relatively small 
scenario of 100 cells and 300 TPs, the computation time is already about 200 times larger for 
the MIP algorithm compared to the proposed sMM algorithm. For larger setups with 1000 TPs 
the normalized time to solve the MIP was ~ 237 compared to ~ 0.49 for the sMM algorithm, 
which is an approximate 488 fold reduction in the computation time. We emphasize that the 
simulated scenarios are small and the computation of the MIP solution becomes infeasible in 
practical scenarios. Already for a network with 200 cells and 10,000 TPs, the sMM algorithm 
provided a solution in about 13s, whereas the MIP algorithm could not find a solution within 
one hour. Compared to the cCZ algorithm the proposed sMM algorithm takes longer time due 
to the lower complexity heuristic used in the cCZ algorithm. For a scenario of 300 TPs the 
average computation time is about 22 times larger for the sMM algorithm and with 1000 TPs 
it is about 43 times larger. However, with typical values of less than Is, the computation time 
is still reasonably small to allow for an online implementation. Considering the advantages in 
energy savings, as seen from Fig. [H the proposed sMM algorithm presents a good trade off 
between computation time and energy savings. 

C. Cells with dijfevent sources of energy consumption 

In contrast to other approaches to the problem of energy-efficient network topology control, 
our optimization framework can easily deal with heterogeneous networks in which cells have 
different static and load-dependent energy consumptions in dH). In other words, the proposed 
sMM algorithm can cope with different energy consumption models of cells. It can select those 
network configurations that exhibit as low overall energy consumption as possible. To illustrate 
the impact of different energy consumption models on the optimization result, we start by varying 
the static energy consumption of all cells, while keeping the load-dependent energy consumption 
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number of TPs 

Fig. 2: Comparison of normalized computation time to obtain results with the sMM algorithm, 
the cCz algorithm and the direct solution of the MIP problem. Normalization with respect to the 
empirical average of the MIP's computation time for 100 cells and 100 TPs over 100 realizations 


fixed. Later in this section, we show the impact of the load-dependent energy consumption by 
changing the weight of the load dependent part relative to the static part. 

To study the impact of the static energy consumption of cells e*, in the following simulations 
we use single cell omnidirectional base stations, and we set the load-dependent part for all 
cells and the common static part at base stations to zero f{pi) = 0 and q = 0. The static 
energy consumption of half of the cells is varied, while the static energy consumption of the 
other half remains unchanged. We refer to the cells with standard fixed energy consumption as 
type 1, while type 2 is used to refer to cells with a varying energy consumption. The energy 
consumption of type 2 cells is specified relative to that of type 1 cells. More precisely, an energy 
consumption relation of (3 = 0.5 means that if Cj = 780W for type 1 cells, then c* = 390W for 
type 2 cells. The results for a scenario consisting of 100 cells and 100 TPs are shown in Fig. [3l 
The simulation confirms the ability of our optimization framework to incorporate different static 
energy consumptions. When all cells consume the same amount of energy {{3 = 1), the algorithm 
makes no difference between type 1 and type 2 cells. The energy consumption of type 1 and 
type 2 cells is roughly the same indicating that equally many type 1 and type 2 cells are active 
in the obtained solution. In contrast, if type 2 cells consume less energy than type 1 ({3 < 1), 
then the algorithm prefers to deactivate type 1 cells, while attempting to keep type 2 cells active. 
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A - sMM type1 
X - MIP typel 
-B— sMM type 2 
V - MIP type 2 
0 - total sMM 
-S — total MIP 


p = energy consumption type 2 BS / energy consumption type 1 BS 


Fig. 3: Fraction of active type 1 and type 2 cells in the final solution obtained with sMM and 
MIP. Deployment uniformly at random for type 1 and type 2 eells. 


Obviously, if /3 > 1, the situation is reversed in the sense that, if possible, type 2 eells are 
preferably seleeted for deactivation. 

The differentiation beeomes even more evident for eell deployments, where type 1 and type 
2 eells are eo-loeated. In sueh a ease, two eells of different type are loeated at the same site 
and are “exehangeable” with respeet to the serviee provided to the TPs (reeall that we use 
omnidireetional eells in these simulations). In other words, if a TP is assigned to a loeation with 
two eo-loeated eells, then it does not matter whieh eell is used to provide the serviee to the 
TP. This implies that the deeision whether to deaetivate a eell or not should depend only on the 
energy eonsumption of this eell in relation to its eo-loeated eelfl The simulations with sueh a 
deployment are shown in Fig. S where we see that, for (3 < 1, there is no aetive eell of type 
1, while, for > 1, type 2 cells eonsume more energy and the simulations eonfirm that the 
algorithms clearly prefer to activate type 1 cell. 

To obtain insight into the impaet of the load-dependent energy eonsumption, we fix the static 
energy consumption of a single-cell omnidireetional base station to be e* = 780W and c; = OW, 
and we vary the load-dependent energy eonsumption fi{p) = 564c'pj by letting d take values 
on d G {0,1,10}. For an inereasing number of TPs, Figl5] shows the fraetion of aetive eells, 
while the normalized network energy eonsumption is shown in Figjbl First we observe that 
the network energy eonsumption always inereases with an inereasing number of TPs, whieh is 
in faet no surprise. Moreover, the fraetion of aetive eell inereases with d for both the sMM 


Even though such setups are unlikely in practice, we use it for reasons of illustration. 
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p = energy consumption type 2 BS / energy consumption type 1 BS 


Fig. 4: Fraction of type 1 and type 2 cells in the final solution obtained with sMM and MIR 
Deployment of type 1 eells uniformly at random and type 2 eells are eo-loeated with type 1. 



Fig. 5: Fraetion of aetive cells for different dynamie energy eonsumption e’ with inereasing 
number of TPs. 


algorithm and the MIP algorithm. An examination of the objeetive funetion in [5^ shows that this 
is what we expeet beeause if the ratio of the load-dependent energy eonsumption beeomes larger 
relative to the statie one, then the algorithm tends to inerease the fraetion of aetive eells for an 
improved load balaneing in order to keep the load of eaeh aetive eell at a relatively low level. 
In other words, instead of deaetivating as many eells as possible to minimize the statie energy 
eonsumption, the algorithm deaetivates the eells to find the best possible balanee between the 
statie and load-dependent energy eonsumption. This ean be observed in FigjSl where we ean 
see that the higher is the load-dependent energy eonsumption (whieh is refleeted by c' > 0), 
the more eells are aetivated under both the sMM algorithm and MIP algorithm. In partieular, if 
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Fig. 6: Normalized network energy consumption for different dynamic energy consumption 
weights c’ with increasing number of TPs. Normalization with respect to the energy consumption 
when all cells are active. 


d = 10, then the fraction of active cells is significantly increased compared with the situation, 
in which the load-dependent energy consumption is negligible (c' = 0). 

D. Alternating sMM algorithm 

We now study the performance of the alternating sMM algorithm presented in Sec. |Vl The 
standard simulation parameters are used with a total number of Z = 10 iterations. To show the 
effect of different TP requirements, we performed simulations under our standard simulation 
setup for different mean data rates Hd at TPs and, for each mean data rate, we used 100 different 
realizations of the simulation scenario. The initial link spectral efficiency is computed based 
on the worst-case interference according to Eq. (O. Our goal is to show the huge potential for 
energy savings when the actual load is estimated as in Alg. [3l instead of assuming the worst- 
case interference scenario, which corresponds to the full-loaded system (see Definition [5]). The 
outcome of the simulation is depicted in Fig. |71 which includes the 95% confidence level and 
shows the normalized network energy consumption with respect to the energy consumption when 
all cells are active. We can see that the application of the techniques from Sec. |V] leads to a 
significant reduction of the normalized network energy consumption for both the sMM algorithm 
and the optimal MIP solution. Furthermore, the largest reduction was always observed after the 
first iteration, which shows that the worst-case interference assumption is very conservative and 
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Fig. 7: Normalized network energy eonsumption when applying Alg. [3] when using the sMM 
algorithm eompared to the using the MIP solution in eaeh iteration. Normalization with respeet 
to the energy eonsumption when all eells are aetive. 


the load estimation may lead to eonsiderable performanee gains. 

VII. Concluding Remarks 

We have introdueed an optimization framework for enhaneing the energy effieieney of eel- 
lular networks. In wireless systems, problems of this type are hard to solve beeause of their 
eombinatorial nature and the nontrivial interferenee eoupling among eells. Indeed, even with 
a simplifying assumption of the worst-ease interferenee, the energy saving problem is a mixed 
integer programming problem that is strongly related to the bin-paeking problem, whieh in turn is 
known to be NP-hard. As a result, we eannot expeet to find optimal solutions quiekly, we foeused 
in this study on fast sub-optimal heuristies. Unlike many existing approaehes in the literature, 
the proposed methods ean naturally eonsider both the dynamie and statie energy eonsumption 
of base stations with multiple eells in heterogeneous networks. 

In the first proposed heuristie, we relaxed the mixed integer programming problem to a form 
suitable for the applieation of majorization-minimization teehniques. The resulting algorithm 
requires the solution of a series of linear programming problems that ean be effieiently solved 
with standard mathematieal solvers. Therefore, it ean be applied to large-seale problems, and it 
is also suitable for online operation. One limitation of this first method is that it uses the worst- 
ease interferenee seenario, so it ean be too eonservative in terms of energy savings. To address 
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this limitation, we also proposed a two-step alternating approaeh that obtain aeeurate values 
of the speetral efficieney of links by using the framework of standard interferenee functions. 
Simulations show that the proposed fast heuristics are able to obtain network configurations that 
are competitive in terms of energy consumption against optimal algorithms. 
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Appendix 

Here we briefly summarize the majorization-minimization (MM) algorithm [fTTll . which can 
be seen as a generalization of the well- known expectation-maximization (EM) algorithm. The 
presentation that follows is heavily based on that in the study in |[39l (see also [|TI|, [l34ll ). 

Suppose that the objective is to minimize a function : Tf —)■ M, where X C Assume 
that there exists a solution to this optimization problem, and let x* G be a global minimizer 
of h] i.e., h{x*) < h{x) for every x E X. Unless h has a special structure that can be exploited 
(e.g. convexity), finding x* is computationally intractable in general [|4^ . Hence, we typically 
have to content ourselves with generating a sequence of vectors with non-increasing objective 
value. To this end, we can use the majorization-minimization (MM) technique, which drives h 
downhill with the help of a majorizing function : A’ x A" —>■ M. In more detail, we say that g 
is majorizing function for h if it satisfies the following properties: 

C.l g majorizes h at every point in X, i.e. 

h{x.)<g{x,y), Vx,yGA’, (16) 


C.2 g and h coincide at (x, x) so that 

h{x) = 5 f(x,x), Vx G X. 


(17) 
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By starting from a feasible point G X , the MM algorithm generates a sequence C 

X with monotone decreasing function values according to (we assume that the optimiza¬ 

tion problems have a solution) 

g arg min 5 f(x, x^”^) . (18) 

xGA:’ 

Irrespective of the choice of g , we can easily verify monotonicity of the objective value with 
the help of (fT^ . (fTTI) and (fTSl) : > 5 f(x(”+^\ = 

^(x(n-i-i)) xherefore, since the function h is bounded below when restricted to X by assumption, 
we can conclude that h { yS ^^) ^ c G M for some c > /^(x*) as n — )■ cx). However, we emphasize 
that this in general does not imply the convergence of the sequence {x^^^}. 

The choice of the function g is problem dependent, but it should be sufficiently structured 
in order to make the optimization problem in (fT^ tractable. In particular, in our study we deal 
with concave and continuously differentiable functions h . In such cases, a natural choice for g 
satisfying (fT^ and (fTV]) is 

^(x,y) = /z(y) + VMyr(x-y). (19) 

This particular choice is common in, for example, sparse signal recovery ll^ . 

Remark 5: We note that, instead of solving the optimization problem in (fT^ exactly, it is 
sufficient for the monotonicity of the sequence that 5 f(x*^’^+^\x^"^)) < 

for every n G N. This observation is relevant if the right-hand side of (fTSl) can only be solved 
asymptotically, in which case the iteration can be truncated whenever the above inequality is 
satisfied. 
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